Effective viscosity of microswimmer suspensions 
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The measurement of a quantitative and macroscopic parameter to estimate the global motility of 
a large population of swimming biological cells is a challenge Experiments on the rheology of active 
suspensions have been performed. Effective viscosity of sheared suspensions of live unicellular motile 
micro-algae (Chlamydomonas Reinhardtii) is far greater than for suspensions containing the same 
volume fraction of dead cells and suspensions show shear thinning behaviour. We relate these 
macroscopic measurements to the orientation of individual swimming cells under flow and discuss 
our results in the light of several existing models. 

PACS numbers: 47.63.-b, 47.57.-s, 47.50.-d 



In nature, organisms that can propel themselves in a 
fluid medium are ubiquitous. While larger organisms, 
such as fish, use inertia in their motion, micro-organisms 
like spermatozoa, micro-algae or bacteria, move at low 
Reynolds number, where viscous forces dominate over the 
effects of inertia. A recent and currently unresolved is- 
sue involves understanding the hydrodynamics associated 
with the individual or collective motion of microswim- 
mers through their fluid- mediated interactions Mi- 
croswimmer suspensions have been shown to lead to com- 
plex dynamics such as the so-called weak turbulence or 
bioconvection phenomenon [l-Q. Such active suspen- 
sions are made of self-propelled particles that create a 
force multipole either at the front of the body, in which 
case they are called pullers, or at the back, in which case, 
they are pushers. The flow induced by the force multi- 
pole is responsible for hydrodynamic interactions which 
are expected to have a dramatic effect on dynamics and 
particularly on the rheology of the suspension. 

Recently, theoretical efforts have been made to model 
the effective viscosity of active suspensions. Stokesian dy- 
namic simulations by Ishikawa and Pedley [5] show a dif- 
ference in effective viscosity for suspensions of swimming 
spherical 'squirmers' in a gravity field. Haines et al. [6| 
analytically showed that swimming results in a change 
in viscosity if the orientation distribution of swimmers is 
assumed to be anisotropic. More recently, Saintillan [?) 
used a simple kinetic model to study the rheology of a di- 
lute suspension of self-propelled particles in a shear flow. 
He showed that suspensions of pullers exhibit increased 
effective viscosity compared to passive suspensions, while 
pusher suspensions exhibit a significant decrease in vis- 
cosity due to the motile activity; these results are consis- 
tent with previous predictions of Hatwalne et al. [sl who 
used coarse-grained active hydrodynamics to predict the 
rheology of active suspensions. Depending on swimmer 
types (puller or pusher), cell shape (spherical or ellip- 
soidal) and locomotion mechanisms, models can lead to 
very diverse results in terms of rheology. In order to bet- 
ter understand the effects of motility on viscosity, exper- 
imental work addressing the rheology of microswimmer 
suspensions is necessary. 



Very recently, Sokolov and Aranson ^ measured the 
microrheology of suspensions of pusher-like bacteria. 
They found, as predicted by f?!, (sj, that the effective 
viscosity of such active suspensions decreases in com- 
parison to passive particles at the same volume fraction. 
To our knowledge, no previous experimental data exist 
for puller-type micro-swimmer suspensions for which an 
increase of effective viscosity is predicted 0, This 
paper presents the first direct experimental macroscopic 
measurement of the effective viscosity of puller type mi- 
croswimmer suspensions: Chlamydomonas Reinhardtii 
(CR), a 10 fiia motile unicellular alga. Rheological mea- 
surements performed on this system show a clear in- 
crease in effective viscosity when compared to a dead 
cell suspension. Shear thinning behaviour is also mea- 
sured. Based on models elaborated for dilute suspen- 
sions, we provide an interpretation of our results, invok- 
ing the anisotropic distribution of live cell orientations 
within the shear flow. We then discuss two hypotheses 
that could account for the origin of this anisotropy: a 
gravity torque or an effective elongated aspect ratio due 
to flagella beating. 

CR micro-algae [10] is a genus of green alga. It is a 
bi-flagellated unicellular organism. Chlamydomonas is 
used as a model organism for molecular biology, espe- 
cially for studies of flagella motion, chloroplast dynamics, 
biogenesis and genetics. They are spheroidal in shape 
with two anterior flagella. Their back-and-forth move- 
ment produces a jerky breast stroke with a mean speed 
oiV 4:0fim/s in a water-like viscous medium. Since the 
cell radius is R ^ bfim^ Brownian motion is negligible. 
The swimming direction of the cells can be directed by 
stimulus gradients: a phenomenon known as taxis, such 
as chemotaxis, rheotaxis or phototaxis. Gradients are not 
used in our experiments in order to avoid any external 
tropism on the motility. Wild -typ e strains were obtained 



from the IBPC lab in Paris [ll|. Synchronous cultures 
of CR were grown in a Tris- Acetate Phosphate medium 
(TAP) using a 14/10 hr light/dark cycle at 25^C Cul- 
tures were typically grown for two days under fluores- 
cent lighting before cells were harvested for experiments. 
These cultures were concentrated up to a typical volume 
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FIG. 1: a: Two dimensional trajectories of 50 swimming cells. Trajectory duration is 20 seconds. Trajectories starting positions 
have been all shifted to the origin, b: Measured mean square displacement (MSD) of cells as a function of time. Averages have 
been extracted from 9Hz sequences of about 700 cells trajectories. MSD fits well to a persistent random walk process with a 
correlation time tv = 3.5 ± 0.1 s and a coefficient D — 995 =b 20 fim'^ .s~^{see text). Inset: Log-Log plot, c: Measured mean 
velocity V of cells as a function of the inversed viscosity of the medium. The line represents a linear fit. Inset: Viscous drag 
force should be proportional to the product of the medium viscosity 770, the mean radius R and the velocity V of the cells. This 
product rjoRV is shown to be independent of the medium viscosity. Error bars are smaller than the symbols. 



fraction of 20 to 30% by centrifuging for 20 minutes at 
900^ and re-suspending them in a fresh culture medium 
to achieve the volume fractions required. Volume frac- 
tions were measured using hematocrite capillaries and a 
Neubauer counting chamber. 

We started by characterizing cell motility. Microscopy 
imaging of chlamydomonas suspensions was carried out 
on an Olympus inverted microscope coupled with a Sen- 
sicam camera used at frame rates up to 30 Hz. Cham- 
bers made of glass and 240 /im spacers were coated with 
bovine serum albumine in order to reduce cell adhesion. 
Cells were imaged in bright field using xlO magnifica- 
tion. Two dimension trajectories of cells were recorded, 
typically for a few tens of seconds. Particle tracking 
was performed using IDL The cell trajectories ob- 
served are correctly modelled by a persistent random 
walk [13I, [31. Short term correlations in the direction 
of movement are observed: the swimmer then presents 
an almost fixed direction during a characteristic time ty. 
This short time stage corresponds to a ballistic regime 
characterized by a velocity V. The ballistic regime ends 
when swimmers make a turn. Note that the reorien- 
tation process has been studied recently by R. Gold- 
stein [15| and is due to a dephasing of the two flagella 
for a short time. Hence the mean square displacement of 
cells is well described by a random walk behaviour with 
a persistence length £p of approximately ten cell diame- 
ters: < + ^2 ^2 ^/^^ _ Q 5 ^2 _ [-2t/ty)] 

where <> represents an ensemble average over more 
than a thousand independent measurements [12]. For 
t « ty, < ^ y'^ >r^ V'^t^ with the ballistic veloc- 
ity V = ip/ty and for longer times a random walk is 
observed, < + >~ 4Dt with D = ll/Uy. 

The algae velocity dependency on viscous drag was 
measured by adding a small amount of short chain dex- 
tran {Mw = 20000) to the culture medium. This al- 



lowed variation of the viscosity of the medium between 
1 and 2.5mPa.s. We checked that swimming velocity is 
inversely proportional to bath viscosity 770 (fig. [Ij). This 
suggests that the stall force (the force needed to stop 
the swimming cell), which is proportional to the product 
rjoRV^ remains constant (fig. [He inset) 

Let us now turn to the actual viscosity measurements. 
The effective viscosity of algae suspensions is measured 
on a Bohlin Gemini 150 rheometer equipped with cone- 
plate geometry (cone angle = 2^, diameter = 60mm). 
Steady shear measurements were made at T = 20^C If 
a shear stress a is imposed and the associated shear rate 
7 measured, effective viscosity is r^eff = cr/7. Samples 
of different volume fractions were prepared and the ef- 
fective viscosity measured for both live and dead cells. 
The cells were checked for motility after the rheometric 
measurement. 

The effective viscosity r^eff of active suspensions is 
found to decrease with the shear rate (fig. [2la). As dis- 
cussed below, it reveals competition between the shear 
rate and cells orientation in the flow. Note, that viscosity 
does not reach a plateau value. This suggests that other 
effects like segregation might play a role at higher shear 
rates. However, a clear global decrease of rjeff is observed 
from the maximum measured value at 7 ^ 4^"^. We then 
investigated the dependence of effective viscosity on the 
volume fraction of the suspension. To do so, viscosity 
was measured at a given shear rate of 5s~^, which is suf- 
ficiently high for rheometer resolution but low enough to 
allow viscosity to be affected by motility. 

Fig. Ob shows the relative viscosity (77eff — ^o)/^o of 
live and dead cell suspensions as a function of the volume 
fraction. In both cases, viscosity is an increasing function 
of the volume fraction as it is for passive beads. Remark- 
ably, the effective viscosity of swimming cell suspensions 
is quantitatively larger than the viscosity of dead cell sus- 
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FIG. 2: a. Effective viscosity of chlamydomonas suspensions as a function of shear rate. Data are shown for different volume 
fractions of the suspension and star symbols represent the viscosity 770 of the culture medium. The lines are shown for ease 
of viewing only. b. Relative viscosity of microswimmer suspensions (measured at shear rate = 5s~^) as a function of volume 
fraction. Solid symbols represent live cell data and crossed symbols represent dead cell data. Measurements are fitted to 
equation [1] using (/)m = 0.62 and a = 2.5 (dead cells) and 4.5 (swimming cells). 



pensions (up to a factor of 2 for a 15% volume fraction). 

The effective viscosity r^eff of a suspension of passive 
spherical particles in a solvent of viscosity 770, depends 
on its volume fraction cf). Krieger and Dougherty's semi- 
empirical law [IgI is shown to provide a reliable descrip- 
tion of the measurements: 

r?efr = 770(1-;^)-"^-. (1) 

(Pm 

Here, (j)m is the maximal packing volume fraction. 
For a dilute regime, where (j) <C 1, eq.([T|) reduces to 
Veff ^ ^0(1 + where a is known as Einstein's in- 
trinsic viscosity {a = 2.5 for passive and spherical parti- 
cles) [13. 

To quantify the effect of motility on the effective viscos- 
ity of cell suspensions, we fitted the measurements using 
eq. ([!]). We set the maximal packing volume fraction to 
0.62 and left a as a free parameter. This resulted in a 
large difference between active and passive suspensions. 
We found a = 2.5±0.1 in the case of dead cells, which 
clearly behave as passive particles. However, in the case 
of swimming cell suspensions, we obtained a = 4.5 ± 0.2. 
Swimming cells induce a quantitative increase in the ef- 
fective viscosity of the suspension. 

In order to understand this increase in effective viscos- 
ity, we conducted complementary experiments which con- 
sisted of imaging the cells while subjected to a shear flow. 
Dead cells show a regular rotation at an angular velocity 
close to 7/2, as would passive spherical particles. We ob- 
served that swimming cells behave very differently: they 
resist the flow rotation for most of the time and even- 
tually flip very rapidly. Fig. [3^ and b show picture se- 
quences extracted from a fast-image film. High frequency 
acquisition {500Hz) allows us to determine whether or 
not an alga is swimming by looking at the beating of 
flagella. In the case of a dead cell, fiagella only move 
because of thermal agitation, whereas flagella of swim- 
ming cells beat at about bOHz. Fig. [3] shows 50i^2; time 
sequences of cells subjected to a 10s~^ shear rate. The 



swimming cell spends about two thirds of the period in 
the xO^— plane and flips during one third of the period 
whereas the dead cell rotates at a constant velocity close 
to 7/2. 

Discussion. A puzzling question arises from our experi- 
ments: why should live cells and dead cells with the same 
aspect ratio and the same volume fraction respond so dif- 
ferently in a flow? In order to answer such a question, we 
propose some hypotheses based on existing models which 
might explain the increase of effective viscosity according 
to individual dynamics of live cells. As shown on fig.[TJb, 
the typical time ty after which a chlamydomonas changes 
its direction of motion is about 3 to 4 seconds. This is 
much longer than the typical time associated with the 
shear flow: tf = 7"^ ^ 0.2s. This means that live algae 
resist the flow rotation (fig. 3- A) during approximately 
70% of the time ty by being aligned in the flow. The 
fast flip between two aligned positions is indeed very fast 
in comparison to the time the cell spends parallel to the 
flow. After a time ty, a reorientation occurs but the live 
alga is then aligned in the flow once again and so on. 
On the contrary, a dead cell rotates at a uniform angular 
velocity 7/2 like a passive particle. 

As proposed by Ramaswamy [8] and Saintillan 0, 
when active particle distribution is anisotropic (i.e. 
ahgned in the flow) it gives rise to a modified effective 
viscosity which is enhanced in the case of puller type cells 
(and reduced in the case of pusher type cells [9|). Fol- 
lowing 0, m a swimming stress tensor can be defined: 
T,^ = cro[< p.p > —1/3]- ctq is the amplitude of the dipo- 
lar forces exerted by the ceh to swim (for pullers ctq > 0). 
The vector p represents the orientation of the cell and in 
our case, it is oriented from the centre of each CR to the 
middle of the flagella anchor points. Active suspensions 
with Q =< p.p > —1/3 7^ are thus equivalent to an 
orientational ordered suspension with a non-zero devia- 
tronic swimming stress contribution. This contribution 
to the effective shear viscosity is such that r]eff = ^xy/i- 



4 



Since for puller ctq > 0, effective viscosity is increased. 
These models although restricted to dilute regimes allow 
us to relate the observed cell alignment in flow to the 
resulting increase in viscosity. 

We can think of two different hypotheses to explain 
the individual dynamics observed (figlSlb): either grav- 
ity effects or large effective aspect ratio due to flagella 
beating. Gravity can exert a torque on cells since the 
centre of mass of CR is known not to coincide with its 
geometrical centre Q. Ishikawa and Pedley [HI showed 
that a suspension of bottom heavy squirmers results in 
an increased effective viscosity. In 1970, Brenner [l8| an- 
alytically predicted the effect of gravity on the rheology 
of a dilute suspension of inhomogeneous particles. He 
obtained for non rotating particles an intrinsic viscosity 
a = f + f . This value is very close to our measurement 
of a. Since dead cells are observed to tumble like passive 
spherical particles, this would suggest that mass distribu- 
tion is homogeneous unlike live cells. However, we have 
no experimental evidence that this is the case. If the 
shear rate 7 is increased, the constant torque To exerted 
by gravity will eventually become smaller than the torque 
necessary to stop cell rotation Tf = 4:7Tr]R^'y. This would 
explain the shear thinning behaviour observed and also 
predicted by Brenner [Hj in figlSla. 

Another hypothesis that would explain the alignment 
of live cells is that their systematic flagella beating results 
in an effective hydrodynamic aspect ratio larger than that 
of their quasi-spherical shape [l9| . The observed individ- 
ual dynamics (figlS]) would then correspond to Jeffery's 
orbits of ellipsoidal prolate particles [20] that result in 
an increased effective viscosity. This is also predicted in 
0, [sl where the puller suspension is close to a suspension 
of passive Brownian rods for which it is well known that 
a shear thinning behaviour can be observed and [2l|. 
Note that Brownian contribution is replaced here by hy- 
drodynamic diffusivity due to hydrodynamic interactions 
between swimming CR. 

In this paper, we have experimentally shown that ac- 
tive suspensions of puller type microswimmers present 
a dramatic increase in effective viscosity. We have cor- 
related these macroscopic rheological measurements to 
the individual dynamics of cells under shear flow: swim- 
ming cells tend to resist the flow rotation unlike dead 
cells. More experiments are needed to explain the origin 
of this phenomenon: a gravity torque exerted on the cell 



or an effective large aspect ratio due to flagella beating. 

We thank C. Nizak for his help in the choice of the 
system and C. Verdier for the rheometry. This work is 
financed by the Rhone- Alpes Region through the Cible 
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FIG. 3: A. Dead cell in a 10 shear flow experiencing tumbling around the x— axis. B. Live chlamydomonas swimming 
the same shear flow. Arrows indicate cells flips and time between pictures is 20ms. The bar represents 10 micrometers. 
Schematic view of the flow cell, where Oz is the direction of observation, yOz is the shear plane. 



